Reheating and leptogenesis in a SUGRA inspired brane inflation 
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We have studied extensively phenomenological implications in a specific model of brane inflation 
driven by background supergravity [H, via thermal history of the universe and leptogenesis per- 
taining to the particle physics phenomenology of the early universe. Using the one loop corrected 
inflationary potential we have investigated for the analytical expression as well as the numerical 
estimation for brane reheating temperature for standard model particles. This results in some novel 
features of reheating from this type of inflation which have serious implications in the production 
of heavy Majorana neutrinos needed for leptogenesis through the reheating temperature. We have 
also derived the expressions for the gravitino abundance during reheating and radiation dominated 
era. We have further estimated different parameters at the epoch of phase transition and revealed 
their salient features. At the end we have explicitly given an estimate of the amount of CP violation 
through the effective CP phase which is related to baryon asymmetry as well as gravitino dark 
matter abundance. 



I. INTRODUCTION 



It is now well accepted that the post big bang universe 
passed through different phases having two-fold signif- 
icance - phenomenological and cosmological. One of the 
significant phases, namely, reheating 0] plays the pivotal 
role in explaining production of different particles from 
inflaton/ vacuum energy. As we look back in time reheat- 
ing was completed within the first second (and probably 
much earlier) after the big bang. At that time nucleosyn- 
thesis [4|, or the formation of light nuclei occurred. Parti- 
cle physicists as well as cosmologists have a clear picture 
of this hot big bang phase because ordinary matter and 
radiation were driving it and also the physical processes 
that characterize it involve terrestrial physics. On the 
other hand the mysterious force that drives the inflation- 
ary phase is conventionally described by a scalar field, 
named inflaton which oscillates near the minimum of its 
effective potential and produces elementary particles Q . 
These particles interact with each other and eventually 
they come to a state of thermal equilibrium at some ar- 
bitrary temperature T. This process completes when all 
the energy of the classical scalar field transfer to the ther- 
mal energy of elementary particles. Since long theoretical 
physicists have been investigating reheating as a pertur- 
bative phase @, or one in which single inflaton quanta 
decayed individually into ordinary matter and radiation. 
The recent theoretical studies have shown that in many 
cases the decay occurs through a non-perturbative pro- 
cess , in which the particles behave in an ordered man- 
ner. Non-perturbative processes involved at reheating 
are extremely more efficient than the perturbative ones 
Q and often more difficult to investigate in practice. In 
short there is no existence of a complete theory which ex- 
plains non-perturbative effects during reheating for the 
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total time scale. 

Besides production of gravitinos during perturbative 
reheating [SHIll its decay plays a significant role in the 
context of leptogenesis pjl, Qjl ■ More precisely two types 
of gravitinos are produced in this epoch - stable [lij and 
unstable [Hj]. Stable ones and decay products of unsta- 
ble ones directly or indirectly stimulate the light element 
abundances during big bang nucleosynthesis. Most im- 
portantly the unstable one has important cosmological 
consequences out of which the major one directly affects 
the expansion rate of the universe (l6| . In order to ex- 
plain cosmological consequences at a time by a single 
physical entity, it is customary to explain everything in 
terms of gravitino energy density which is directly pro- 
portional to the gravitino number density or gravitino 
abundance. This gravitino abundance is obtained by 
considering gravitino production in the radiation dom- 
inated era following reheating [l7| . Gravitinos are orig- 
inated through thermal scattering [ll], [HI in the early 
universe and are usually related to the reheating temper- 
ature (T reh ). Particle physics phenomenology usually 
requires that under instantaneous decay approximation 
19] reheating temperature (T reh ) is maximum during re- 
heating. 

In the present article we have studied extensively re- 
heating phenomenology and leptogenesis in a typical 
brane inflation model which was proposed earlier by us 
[lj . Precisely, the model includes one loop radiative cor- 
rection in the framework of local brane version of the 
supersymmetric theory i.e. N = 1, D = 4 SUGRA which 
is derived from the background N = 2, D = 5 SUGRA 
in the bulk (for details please refer to [l|). In the present 
article, our prime objective is to investigate for both the 
analytical and numerical expression for brane reheating 
temperature in high energy limit for standard model par- 
ticles followed by a detailed investigation for gravitino 
phenomenology and leptogenesis. Here, and throughout 
the rest of the article, by high energy limit implies that 
the total energy density is very high with respect to the 
brane tension as mentioned in our earlier paper Q. As it 
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will be revealed, the scenario is somewhat different in the 
context of reheating from brane inflation which results in 
novel features worth studying in details. This has seri- 
ous implication for the production of the heavy Majorana 
neutrinos needed for leptogenesis We further esti- 

mate different parameters related to reheating and lep- 
togenesis at the epoch of phase transition [20] . Last but 
not the least we have given an estimate of CP violation 
which is the indirect evidence of the baryon asymmetry 
and connected with gravitino dark matter abundance. 



II. REHEATING PHENOMENOLOGY ON THE 
BRANE FOR SU(2) L <g) U(1) Y 

A. Model Building from background supergravity 

For systematic development of the formalism, let us 
briefly review from our previous paper [1| how one can 
construct the effective 4D inflationary potential of our 
consideration starting from N = 2, D = 5 SUGRA in the 
bulk which leads to an effective N = 1, D = 4 SUGRA in 
the brane. Considering the fifth dimension is compact- 
ified on the orbifold S 1 /Z2 of comoving radius R, the 
N = 2, D = 5 bulk SUGRA is described by the following 
action 



and e4, 5D Planck mass M5 and the brane tension A and 
two constants /3 and 60 comes from the metric structure. 
Here C(T, T^) represents an arbitrary function of T and 
Tt. This leads to an effective N = 1,D = 4 SUGRA in 
the brane with the F-term potential 



V = Vf = exp 1 



( Kj&£) \ \( dW 1 f dK \ w\ 

\ TP J \d^ a ~ r J TP J 



K 



raw 1 / dK 



w 

TP 



,\W\ 2 
~KP~ 



(4) 

Here \E' a is the chiral superfield and <fi a be the 4D com- 
plex scalar field. In this context the Kahler poten- 
tial is dominated by the leading order term i.e. K = 
^2 a 4>a4' a - The superpotential in eqn(0} is given by 
W = T,n=o D nW n (<p a ) with the constraint D = 1. Ex- 
panding the slowly varying inflaton potential around the 
value of the inflaton field along with Z2 symmetry the 
required renormalizable one-loop corrected inflaton po- 
tential turns out to be 
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Here the sum includes the walls at the orbifold points 
Ui = (0,nR) and 5-dimensional coordinates x m = (x a ,y), 
where y parameterizes the extra dimension compactified 
on the closed interval [— wR, +ttR}. Written explicitly, 
the contribution from bulk SUGRA in the action 
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and including the radion fields (x, T, T') the ef- 
fective brane SUGRA counterpart turns out to be 
5(y)L 4 = -e (5) A(y)[(d a 0(d a <f>) + ix<t a D a x]. The 
Chcrn-Simons terms can be gauged away assuming cubic 
constraints and Z2 symmetry. Further, S 1 /Z2 orbifold 
setting allows us to express the 4-dimensional part of the 
action (after dimensional reduction) as, 
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27r2M 4 and Da = Ca — ^jjp wnere C4 i s 
negative constant appearing at the tree level. Here A 
represents the energy scale of brane inflation which can 
be expressed in terms of the slow roll parameter rj v ex- 
plicitly derived in 0. For our model A ~ 0.2 x 10 16 GeU 
for the window -0.70 < D 4 < -0.60. 

With this brief review of the construction of the po- 
tential we are now in a position to investigate for its 
phenomenological significances. 

From the knowledge of particle physics it is known that 
during the epoch of reheatin g in flatons decay into diffcr- 
/mn ent P art i c ^ e constituents pi. l2l| are directly related to 
the trilinear coupling of the inflaton field. There might 
be a possibility of collision originated through quartic 
coupling and driven by background scalar field. For ex- 
ample here the contribution from the heavy Majorana 
neutrino comes from the seesaw Lagrangian CMajo — 
—IiIlRiP — ^Ai^ip + h.c, where II and T-L are the lep- 
ton and the Higgs doublets, respectively, and M. is the 
lepton-number-violating mass term of the right-handed 
neutrino. Now using the assumption ^> m a ,m$ ^> 
the total inflaton decay width for the positively and 
negatively charged 4>((f> + , (f>~) scalar fields as well as the 
fermionic field tp (Example: For the heavy Majorana neu- 
trinos the decay process ip ~ * IlH, ip — > IlR predomi- 
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Q = in^R 2 an< ^ * ne 4D Planck 



where P = 
mass Mpl 

1.22 x 10 19 GeU. Here we have introduced the reduced 4D 
Planck mass M = 2.43 x 10 18 GeU, 5D and 4D charge e 5 



where the coupling strength C 
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j-p , and the background scalar field is a. 

Now to construct the thermodynamical observable the 
effective number of particles incorporating relativistic de- 
grees of freedom is defined [H[ as N* = Ng + iN F , 



3 



where N 



B 



J2j N£j. Here N* B repre- 
sents the number of bosonic degrees of freedom with mass 
77i0 <C T and N F represents number of fermionic degrees 
of freedom with mass <C T. Here 'i' and 'j' stand for 
different bosonic and fermionic s pec ies respectively. For 
the phenomenological estimation [23J N* ~ 10 2 — 10 4 and 
for realistic models N* ~ 10 2 — 10 3 . For convenience let 
us express reheating temperature on the brane as 
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where H(T br ) and p(t reh ) be the Hubble parameter and 
energy density during reheating respectively. It is worth 
mentioning that the brane reheating temperature does 
not depend on the initial value of the inflaton field and 
is solely determined by the elementary particle theory of 
the early universe. 



B. Phase transition in brane inflation 

Phase transition in braneworld scenario is weakly first 
order in nature (24| . So it is convenient to write the brane 
reheating temperature in terms of the critical parameters. 
To serve this purpose the critical density and the critical 
temperature or transition temperature can be written as 



p(t c ) = 2\ = 
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which makes a bridge between the phenomenology and 
observation. In the high energy limit 5D Planck mass 
(M5) can be expressed in terms of our model parame- 
ters as M5 = ^/ 64007r z M^" 4 + 4i ? iL(^ Here A 2 repre- 
sents the amplitude of the scalar perturbation defined as 

Most importantly here 
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the subscript * represents here the epoch of horizon cross- 
ing (k — aH) and a represents a dimensionlcss model 
parameter defined as a = 

The major thermodynamic quantities - critical density 
(p c ), critical pressure (P c ), critical entropy (S c ) - and 
the Hubble parameter at the critical temperature (H c ) 
related to the phase transition designated by a four tuple 
critical characteristic set U (07) by the following fashion 
for our model: 



U{crj) : [{ Pc ,P c ,S c ,H c } 
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where we have defined a dimensionless characteristic 
quantity Ayq)*) 



a 4 M 2 



at the horizon cross- 
ing in this context. The above mentioned physical quan- 
tities are function of the critical or transition temperature 



which is defined as 
T r := 



288000 
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with gauge group J := SU(2)p ® U(l)y and the species 
index 7 = 1(B Boson), 2(F =4> Fermion),3(M 
Mixture). 




FIG. 1: Here we have plotted the variation of the critical den- 
sity with respect to the effective inflaton field (f>+ after horizon 
crossing in the domain —0.70 < D4 < —0.60, which explic- 
itly shows the direct connection between the phenomenologi- 
cal and observational sector through the 5D Planck mass M5 
given in eqn®. In addition it confirms the existence of phase 
transition in braneworld scenario in high energy limit. Here 
we have used M = 2.43 x 10 18 GeV. Most importantly the 
energy density is smoothly varying function of inflaton field 
in the critical domain. 



C. Brane reheating temperature 

In this context the reheating temperature can be writ- 
ten [25l | as a one to one mapping (0) in parameter space 
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where Q represents collection of all gauge group which 
supports particle theory. But in this context we are con- 
fining ourselves into the Standard Model regime. So to 
construct a fruitful model of reheating in the context of 
Standard Model gauge group, we rewrite all general prin- 
cipal components in terms of physical degrees of freedom 
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in a compact fashion. We consider a one to one high 
energy mapping Q[~f] in a physical space such that 
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it maps the actual brane reheating temperature (Ti"") to 
its high energy value (T brh ) in the Standard Model gauge 
group J := SU(2) L (g> U(1) Y with Z 7 = (5. 



40 8> 
7 ' 3/ 



(600, ^fS, 320) and 7 = 1(5), 2(F), 3(M). Here 

U 7 C/ [ C 7]0 ( 3[7] ^ for wllich J G G- Most impor- 
tantly the superscript 'br' and 'brh' stands for parame- 
ters before and after high energy mapping respectively. 
Here it should be mentioned that the brane reheating 
temperature incorporates all the effects of heavy Majo- 
rana neutrinos as well as the other fermions and bosons 
through the total decay width T to tai- 

The reheating temperature for different species can 
readily be calculated from our model. For a typical 
value of C4 ~ D4 — —0.7 (consistent with [l|), we 
have: for boson T^ h ~ 7.6 x 10 10 GeV, for fermion 
T b F rh ~ 7.8 x 10 10 GeV and for mixture of species T$ h ~ 
6.5 x 10 10 GeV. This is significantly different from GR 
value T reh < fO 6 — 10 7 GeV and is a characteristic fea- 
ture of brane inflation. 



III. GRAVITINO PHENOMENOLOGY ON THE 
BRANE FOR SU(3) C SU(2) L U(1) Y 

Let us now move on to studying how the self interact- 
ing term of our model is directly related to the leptogen- 
esis through the production of thermal gravitinos which 
is a special ingredient for the heavy Majorana neutrinos 
in the leptogenesis. Let us start with a physical situa- 
tion where the inflaton field starts oscillating when the 
inflationary epoch ends at a cosmic time t — t osc — tf. 
Throughout the analysis we have assumed that the uni- 
verse is reheated through the perturbative decay of the 
inflaton field for which the reheating phenomenology in 
brane is described by the Boltzmann equation [22| 
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Here p r and represent the energy density of radiation 
and inflaton respectively and Y ^ is the rate of dissipation 
of the inflaton field energy density. At t = t osc epoch the 
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= A 4 (^p) exp [-IV(i - t osc )] . 
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It is worthwhile to mention here that the inflaton 
field <fi follows an equation of state similar to radia- 
tion rather than matter i.e. = -y-, where = 

, 4" 



A 4 



1 + (A + *4 In (£))(£) 



Now solving 



Fricdmann equation the dynamical character of the scale 
factor can be expressed as 



a(t) = a osc ^ 



2H osc (t ^osc) 



where we use a specific notation a(t osc ) = a osc . 
Plugging ean p^j) and egn flTSj) in eqn (fT2l we get 
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As a whole phenomenological construction of gravitino 
abundance is governed by the above equation. But 
cqn fTT)) is not exactly analytically solvable. So we are 
confining our attention to the high energy limit where 
the Friedmann equation (|13[) can be approximated as 
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whose solution is given by 
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Now using an physically viable assumption t < 1 
the exact solution of the eqntfTTf in the high energy limit 
can be written as 
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Our intention is to find out the extremum temperature 
during reheating epoch which is one of the prime compo- 
nents for the determination of gravitino abundance. In 
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the braneworld scenario this extremum temperature is 
given by 



nbh 
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and it is less than the reheating temperature in brane 
(T brh ). This phenomenon is different from standard GR 
results [10] where we see that the reheating temperature 
shoots up to a maximum value and it gives the upper 
bound of the reheating temperature. But in the present 
context of brane inflation this situation is completely dif- 
ferent i.e. at first temperature falls down to a minimum 
which fixes the lower bound of the reheating tempera- 
ture and rises to a maximum at the end of reheating 
epoch. Using eqn flTTl) . cqn (j2"Tj) and the thermodynamic 
background of energy density of radiation we can express 
the scale factor in terms of temperature as 



a(T) = 



1-32 



n 32 



if t — tosc — if 

if tosc(— tf) < t < t re h- 



It is worth mentioning that if we break the time scale 
into two parts t osc < t < t ex and t ex < t < t re h, as 
done in GR the scale factor and hence the remaining 
results have same expressions in these two different zones. 
This is in sharp contrast with standard GR results except 
at t = tf, where they have different values in the two 
different regimes. 

Let us now use this phenomenological background to 
derive the expression of the gravitino production during 
two thermal epochs - reheating and radiation dominated 
era. It is well known that gravitinos are produced by 
the scattering of the inflaton decay products 26]. The 
master equation of gravitino phenomeno logy as obtained 
from 'Boltzmann equation.' is given by fl8l |27| 
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where n — is the number density of scatter- 

ers(bosons in thermal bath) with £(3)=1. 20206.... Here 
^totai is the total scattering cross section for thermal 
gravitino production, v is the relative velocity of the in- 
coming particles with (v) = 1 where (...) represents the 



thermal average. The factor 



(Es) 
2 



represents the aver- 



aged Lorentz factor which comes from the decay of grav- 
itinos can be neglected due to weak interaction. For the 
gauge group E := SU(3) C ® SU(2) L (g) U(1) Y the ther- 
mal gravitino production rate is given by, 
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FIG. 2: In the above figure we have plotted the variation of 
the Hubble parameter with respect to dimensionless param- 
eter j n the domain —0.70 < D4 < —0.60, which shows 
the smooth behavior of Hubble parameter except x < i.e. 
a t -fih ^ "J^f ' M° s t importantly here the equality corre- 
sponds to the end of reheating epoch and the beginning of 
radiation dominated era which is the direct outcome of the 
first expression at t ~ tf for the scale factor (a(T)) stated in 
eqn (|22p . The rest of the the part follows the second expres- 
sion given in eqn (|22fl in the interval tf < t < t re h plotted in 
the above figure. Additionally the vertical scale corresponds 
to M = 2.43 x 10 18 GeF. 



where i = 1,2,3 stands for the three gauge groups 
U(1) Y ,SU(2) L and SU(3) C respectively. Here M, repre- 
sent gaugino mass parameters and gi(T) represents gaug- 
ino coupling constant at finite temperature (from MSSM 
RGE) 



i(X) 



(25) 



with bi = 11, 62 = 1, &3 = —3. Here Ci and Ki represents 
the constant associated with the gauge groups with C\ = 
11, C 2 = 27, C 3 = 72 and K x = 1.266,^2 = 1.312, K 3 = 
1.271. 

For convenience let us recast eqn ([23| as 
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where a boundary condition T = T e " 1 ,T = is intro- 



duced. In terms of a dimensionless variable 
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can (|2"6"]) can be expressed as 
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where di 



exact solution of the eqn ([2"51) is given by 



n d (x) = |#v / ^TT(-2 2 F 1 [£;l-di;§;a; + l] 
+ 2 Fi [i;2-di;|;a;+l] +2^1 [§; -dr, f; a: + 1]) , 

(29) 

Using the properties of Gaussian hypergeometric func- 
tion for x >> 1 eqn (|29"]) reduces to the following simpler 
form: 



1 2 1 

) + rFIT rFI) J 
(30) 

Using the boundary condition T — T b £ in eqn (l30l) the 
numerical value of gravitino abundance turns out to be 
nQ(x ex ) — 62.023^3. 
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FIG. 3: In the above diagram we have plotted variation 
of gravitino number density in a physical volume vs scaled 
temperature in braneworlds. Here we have used the funda- 
mental scale (I3 = 4.596 x 10 _44 aM 3 , where a is a dimen- 
sionless number depends on the species of the MSSM gauge 
group E. For an example n — 4 level flat direction content 
QQQL, QuQd, QuLe, uude of MSSM gives a ~ 15.694 in 
the absence of top Yukawa coupling. Most importantly 4D 
effective Planck mass M = 2.43 x 10 18 GeU. From the plot it 
is obvious that the gravitino number density is monotonically 
increasing function of the dimensionless variable except 

ex 

at x < which implies 
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Let us now find out the exact analytical expression for 
the gravitino abundance at reheating temperature T brh 
in the high energy limit. To serve this purpose substitut- 
ing T = T brh in eqn(gni) we get 
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+2-P1 

along with an extra constraint G ( Tbh » 

f( 32 (%r) _1 ) (rff) + rFI) ~~ rFi)}' K is 
convenient to express the abundance of any species 'cr' 
as[10] Yfc = ^ where n;, is the number density of the 
species '6' in a physical volume and 's' is the entropy 
density given by s — ^-N*T 3 . Here the master equa- 
tion, for gravitino can be expressed as 

(33) 

the time-temperature rela- 
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Using eqn (|16p and eqn 
tion can be found as: 
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Eliminating T we get the solution of the master equa- 
tion(l33l) in the radiation dominated era as 



Yg{T f ) = Y^(T br ) + Yg- rad (T f ) 



(35) 



where 
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But in eqnp5p the first term on the right-hand side is 
not exactly computable. As mentioned earlier to find 
out exact expression we have used here the high energy 
mapping. 

In the radiation dominated era the dynamical behavior 
of temperature can be mapped as 
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Using this map we finally have 
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G 



( T f) = (§ 



Y b 

G 



C(3) 
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32 



l 

Tf 



rpbh 



1 

q^brh 



(39) 



( T brh \ n ( T brh \ 



The gravitino dark matter abundance and the baryon 
asymmetry is connected through Y~, 



brh 
G 



N* 



where 

T>(< 1) is the dilution factor and the leading contri- 
bution is given by the interference between the tree 
level and the one-loop level decay amplitudes. Here 
the CP-violating parameter is described as [28| &cp — 

r(^i L n)+r(^i L n*) ~ le™ 2 s1u °cp, wnere m v is tne 
heaviest light neutrino mass, v = 17AGeV is the VEV of 
Higgs and Sep is an effective CP phase which parame- 
terize each entries of the CKM matrix. Particularly 8c p 
acts as a probe of flavor structure in supergravity the- 
ories. The complete wash out situation corresponds to 
V = l. 



3.0xl0 4: 



2.5xl0 4: 




5.x HP 



1.x lO -11 
T(mM) 



1.5x10"' 



2.xl(T 



FIG. 4: Here we have plotted the variation of total grav- 
itino abundance vs temperature in the domain —0.70 < D4 < 
—0.60, which clearly shows that gravitino abundance at zero 
temperature shoots up initially to maximum and then be- 
comes constant with respect to temperature during radiation 
dominated era in braneworld scenario. As mentioned earlier 
we have used the fundamental scale j& = 4.596 x 10 _44 a 
, where a is a dimensionless number which depends on the 
species of the MSSM and M = 2.43 x 10 18 GeV. 



Through out all the numerical estimation we have 



taken decay width ~ 2.9 x 10 3 GeV, mass of the in- 
flaton ui,fi ~ 10 13 Gey, final temperature and time at the 
end of reheating T f ~ 10 6 GeV and t f ~ 1.4 x 10 31 GeV 
respectively. For a typical value of G4 ~ D4 = —0.7 
extremum (minimum) temperature during reheating can 
be estimated as ~ 7.0 x 10 10 GeK This clearly 
shows deviation from standard GR phenomenology [10| 
where the extremum (maximum) temperature during re- 
heating T max ~ 1.3 x 10 12 Gey. Similarly for G 4 ~ 
Z?4 = —0.7 the critical temperature for different parti- 
cle species and gravitino abundance at different temper- 
atures obtained from our model are: for boson T c b — 
3.2 x 10 14 GeU, for fermion T cF ~ 3.3 x 10 14 GeV, for 
mixture of species T c m — 2.8 x 10 14 GeV, at reheating 
temperature Y^(T brh ) ~ 8.1 x 10- 34 Gey- 3 d3 and at the 



end of reheating Y 



2.1 x 10- 13 GeU- 3 (i3- We 



b—rad/rri \ 
" G {T f> 

have calculated all the abundances in the fundamental 
unit of d 3 i.e. rf 3 = 6.594a x 10 n GeV 3 , where a is a 
dimensionless characteristic constant originated through 
the thermal gravitino production rate in the context of 
MSSM. Most significantly for different flat direction con- 
tents the phenomenological parameter is different and 
can be calculated from MSSM RGE flow at the one- 
loop level for that flat direction. To obtain a conser- 
vative estimate of gravitino abundance we have taken 
here gaugino masses Mj — > for all gauge subgroups 
within MSSM. For example the fourth level flat directions 
QQQL, QuQd, QuLe, uude give a = 15.694 for a spe- 
cific choice of the f7V(l), SUl(2) and SUc(3) gauge cou- 
plings gi = 0.56, #2 = 0.72 and 33 = 0.85 respectively ob- 
tained from the universal mSUGRA boundary condition 
and consistent with electroweak extrapolation of the so- 
lution of MSSM RGE flow from the energy scale of brane 
inflation A = 0.2 x l0 16 GeV for our model. The linear de- 
pendence on T hrh makes simple to revise the constraints 
on T hrh based on the lower limit on the gravitino abun- 
dance - the lower bound on T brh is increased by a factor 
of 1.074. Since T b £ cx T brh , T b £ is not affected much. 
Therefore models of leptogenesis that invoke a small T*£ 
to create heavy Majorana neutrinos are not significantly 
affected. Within 55 < N < 70 and T brh ~ 6.5 x 10 10 Gey 
the entropy density changes. As a consequence the total 
gravitino abundance changes according to fig((3]). It is 
easily seen that P = |, S = consistency relations 
are valid in this context. It is worthwhile to mention here 
that in brane pressure and entropy density of the universe 
falls down to a minimum due to the minimum tempera- 
ture during reheating epoch. However during radiation 
dominated era total entropy density is almost constant 
for both the cases. This clearly shows the deviation from 
standard GR phenomenology. Throughout the analysis 
we have not included the effect of exp[— Y${t — t osc )) in 
the energy density of inflaton p$. One might be con- 
cerned that this will lead to inaccuracies close to 
when most of the gravitinos are produced. However if 
one writes p^ ~ a -4 exp(— T^t) ~ t~ 2 exp(— T^t) for 
t >> t b e x then p'^/prj, = —2/t - Y^. Therefore even till 
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close to tbrh = decreases primarily due to the ex- 

pansion of the universe. Furthermore, near it in- 
creases as T" 1 / 2 ~ y/a ~ ta in brane which is again dif- 
ferent from GR phenomenology where T -1 / 2 ~ y/a ~ t% . 
The thermal leptogenesis in the braneworld can take 
place if the lightest heavy neutrino mass lying in the 
range T brh < A4 < T c . This confirms that the upper 
bound of 5D Planck mass M5 < 10 16 GeV (for our model 
M 5 ~ 7.8 x 10 15 GeU for C 4 ~ D 4 = -0.7), which co- 
incides with the leptogenesis bound implied by the ob- 
served baryon asymmetry. It is important to mention 
here that in the standard cosmology, the thermal lepto- 
genesis in supergravity models is hard to be successful, 
since the reheating temperature after inflation is severely 
constrained to be T reh < 10 6 - 10 7 GeV due to the grav- 
itino problem. However, as pointed out in [25| . the con- 
straint on the reheating temperature is replaced by the 
transition temperature in the brane world cosmology. As 
a result the gravitino problem can be solved even if the 
reheating temperature is much higher. In fact, such in- 
flation models are possible but limited and our model 
is also in that category. Here we are using a prefer- 
able value of the heaviest light neutrino mass from atmo- 
spheric neutrino oscillation data m v ~ 0.05eU and for 
sufficient baryon asymmetry the lightest neutrino mass 
M. ~ 10 10 GeU. For complete washout situation (T> = 1) 
in our model the effective CP phase lying within the win- 
dow 2.704 x 10~ 9 < S C p < 2.784 x 10~ 9 , where S C p is 
measured in degree. Most significantly it indicates that 
the amount of CP violation in braneworld scenario is 
very small and identified with the soft CP phase. Conse- 
quently it has negligibly small contribution to K and B 
physics phenomenology. 

IV. SUMMARY AND OUTLOOK 

In the present article we have studied reheating in 
brane cosmology on the background of supergravity. We 
have exhibited the process of construction of a fruitful 
theory of reheating for an effective 4D inflationary poten- 
tial in N = 1 , D = 4 supergravity in the brane derived 
from N = 2,D = 5 supergravity in the bulk We 
have employed this potential in reheating model building 
by analyzing the reheating temperature in the context 
of brane inflation, followed by analytical and numerical 
estimation of different phenomenological parameters. It 
is worthwhile to mention here that we get a lot of new 
results in the context of braneworld phenomenology com- 
pared to standard GR case. Most importantly we get a 
different numerical value of reheating temperature as well 
the extremum temperature compared to the standard GR 
results. Next using the extremization principle we justify 
that the extremum temperature is the minimum temper- 
ature during reheating which again shows deviation from 



standard GR inspired phenomenology. All these facts are 
reflected in the numerical results of the gravitino abun- 
dance in reheating and radiation dominated era. In the 
context of phase transition we also get different numer- 
ical results for different parameters for standard model 
particle constituents. 

We have further engaged ourselves in investigating for 
the effect of perturbative reheating. To this end we pro- 
pose a theory which reflects the effect of particle produc- 
tion through collision and decay thereby showing a direct 
connection with the thermalization phenomena. To show 
this link more explicitly we put forward both analytical 
and numerical expressions for the gravitino abundance 
in a physical volume in the reheating epoch. Next we 
have found out the gravitino abundance in the radiation 
dominated era. Last but not the least we have expressed 
the total gravitino abundance in a final temperature Tf. 
Most significantly the precision level of all estimated nu- 
merical results is the outcome of the AD effective field 
theory which is analyzes with the arrival of lots of so- 
phisticated techniques. 

Apart from the aforesaid success in estimating phe- 
nomenological parameters there are some added advan- 
tages of our model with reheating in brane which are 
worth mentioning. One of the most significant features 
in the context of braneworld phenomenology is the va- 
lidity of leptogenesis for our model which consequently 
shows the production of heavy Majorana neutrinos in 
the brane. 

In future our aim is to search for the signatures of our 
model for domain wall formation [29] linked to the topo- 
logical defects, 'Q-ball' formation [3fJ] connected with the 
non-topological solitons in braneworld, the role of Lee- 
Wick particles in brane reheating and leptogenesis, pri- 
mordial non-Gaussianity, baryogenesis etc. Last but not 
the least the detailed study of quantum pha se transition 
using Monte Carlo simulation technique [311 ] to minimize 
rapid fluctuation [321 ] or oscillation during measurement 
is also an open issue. We expect to address some of these 
issues in near future. 



Acknowledgments 

SC thanks B. Basu, K. Bhattacharya and A. 
Mukhopadhyay for illuminating discussions and Coun- 
cil of Scientific and Industrial Research, India for finan- 
cial support through Junior Research Fellowship (Grant 
No. 09/093(0132)/2010). SP is supported by Alexander 
von Humboldt Foundation, Germany through the project 
"Cosmology with Branes and Higher Dimensions" and is 
partially supported by the SFB-Tansregio TR33 "The 
Dark Universe" (Deutsche Forschungsgemeinschaft) and 
the European Union 7th network program "Unification 
in the LHC era" (PITN-GA-2009-237920). 



9 



[1] S. Choudhury and S. Pal, larXiv:1102.4206l [hep-th]. 

[2] S. W. Hawking and G. F. R. Ellis, Astrophysical Journal, 
152 (1968) 25; S. W. Hawking and Roger Penrose, Pro- 
ceedings of the Royal Society of London, A 314 (1970) 
529. 

[3] L. Kofman, A. Linde and A. Starobinsky, Phys. Rev. 
Lett. 73 (1994) 3195; L. Kofman, A. Linde and A. 
Starobinsky, Phys. Rev. D 56 (1997) 3258; A. Mazumdar 
and J. Rocher, Phys. Rept. 497 (2011) 85; A. Mazumdar, 
larXiv:1106.5408l [hep-ph]. 

[4] S. Buries, K. M. Nollett and M. S. Turner, Phys. Rev. D 
63 (2001) 063512. 

[5] G. Palma and V. H. Cardenas, Class. Quant. Grav. 
18 (2001) 2233; R. Allahverdi and Manuel Drees, 
|arXiv:hep-ph/0210432| M. Basler and B. Kampfer, In- 
ternational Journal of Modern Physics A 9 (1992) 2033. 

[6] R. Allahverdi, R. Brandenberger, F. Y. C. Racine and A. 
Mazumdar, Annu. Rev. Nucl. Part. Sci. 60 (2010) 27; G. 
Mangano, G. Miele, S. Pastor and M. Peloso, Phys. Rev. 
D 64 (2001) 123509; V. H. Cardenas, Phys. Rev. D 75 
(2007) 083512. 

[7] J. H. Traschen and R. H. Brandenberger, Phys. Rev. 
D 42 (1990) 2491-2504; A. D. Dolgov and D. P. Kir- 
ilova, Sov. J. Nucl. Phys. 51, (1990) 172-177; Y. Shtanov, 
J. H. Traschen and R. H. Brandenberger, Phys. Rev. 
D 51 (1995) 5438-5455; D. Boyanovsky, M. Dattanasio, 
H. J. de Vega, R. Holman and D. S. Lee, Phys. Rev. 
D 52 (1995 ) 6805; D. Boyanovsky , H. J. de Vega and 
R. Holman, |arXiv:hep-ph/"970 1304 J. McDonald, Phys. 
Rev. D 61 (2000) 083513; D. Boyanovsky, M. Dattana- 
sio, H.J. de Vega, R. Holman, D. S. Lee and A. Singh, 
|arXiv:hep ph/9505220| 

[8] R. Allahverdi and A. Mazumdar, Phys. Rev. D 76 (2007) 
103526; D. Boyanovsky, M. Dattanasio, H.J. de Vega, R. 
Holman and D. S. Lee, arXiv:hep-ph/9511361 

[9] E. J. Copeland and O. Seto, Phys. Rev. D 72 (2005) 
023506; V. S. Rychkov and A. Strumia, Phys. Rev. D 75 
(2007) 075011. 

[10] R. Rangarajan and N. Sahu, Phys. Rev. D 79 (2009) 
103534. 

[11] M. Bolz, A. Brandenburg and W. Buchmuller, Nucl. 
Phys. B 606 (2001) 518. 

[12] A. Pilaftsis, J. Phys. Conf. Ser. 171 (2009) 012017; S. 
Davidson, E. Nardi and Y. Nir, Phys. Rept. 466 (2008) 
105; W. Buchmuller, R. D. Peccei and T. Yanagida, 
Ann. Rev. Nucl. Par t. Sci. 55 (2005) 311; Y. Nir, 
|arXiv:hep-ph/0702199| W. Buchmuller, D. Bari, P. M. 
Plumacher, Annals Phys. 315 (2005) 305. 

[13] J. M. Frere, F. S. Ling, M. H. G. Tytgat and V. V. 
Elewyck, Phys. Rev. D 60 (1999) 016005. 

[14] W. Buchmuller, K. Hamaguchi, M. Ibe and T. T. 
Yanagid, Phys. Lett. B 643 (2006) 124; T. Asaka, K. 
Hamaguchi and K. Suzuki, Phys. Lett. B 490 (2000) 136- 



146; T. Moroi, H. Murayama and M. Yamaguchi, Phys. 

Lett. B 303 (1993) 289. 
[15] K. Kohri, T. Moroi and A. Yotsuyanagi, Phys. Rev. D 

73 (2006) 123511; L. Covi, M. Grefe, A. Ibarra and D. 

Tran, JC AP 0901 (2009) 029. 
[16] T. Moroi, |arXiv:hep-ph/9503210| 

[17] K. Choi, K. Hwang and H. B. Kim and T. Lee, Phys. 
Lett. B 467 (1999) 211; R. Rangarajan and N. Sahu, 
Mod. Phys. Lett. A 23 (2008) 427; R. Allahverdi and A. 
Mazumdar, JCAP 0610 (2006) 008; L. Boubekeur, K. Y. 
Choi, R. R. de Austri and O. Vives, JCAP 1004 (2010) 
005. 

[18] J. Pradler and F. D. Steffen, Phys. Rev. D 75 (2007) 
023509. 

[19] G. F. Giudice, E. W. Kolb and A. Riotto, Phys. Rev. 
D 64 (2001) 023508; C. Pallis, Astroparticle Physics 21 
(2004) 689-702. 

[20] P. Creminelli, S. Dubovsky, A. Nicolis, L. Senatore and 
M. Zaldarriaga, JHEP 0809 (2008) 036; T. W. B. Kibble, 
Physics Reports 67 (1980) 183. 

[21] D. J. H. Chung, E. W. Kolb and A. Riotto, Phys. Rev. 
D 60 (1999) 063504; L . F . Abbott, E. Farhi and M. B. 
Wise, Phys. Lett. 117 B (1982) 29; J. H. Brodie and D. 
A. Easson, JCAP 0312 (2003) 004. 

[22] E. Kolb and M. Turner , The Early Universe , Addison 
Wesley (1990). 

[23] A. Linde, Contemp. Concepts Phys. 5 (2005) 1. 

[24] J.V. Narlikar and T. Padmanabhan, Annu. Rev. Astron. 
Astrophys. 29(1991) 325. 

[25] R. G. Felipe, Phys. Lett. B 618 (2005) 7. 

[26] A. Ferrantelli, University of Helsinki, Finland, Ph.D. the- 
sis (2010) larXiv:1002.2835l [hep-ph] . 

[27] J. Pradler and F. D. Steffen, Phys. Lett. B 648 (2007) 
224. 

[28] N. Okada and O. Seto, Phys. Rev. D 73 (2006) 063505. 

[29] Y. I. Takamizu and K. I. Maeda, Phys. Rev. D 70 (2004) 
123514; D. Sez and V. J. Ballester, Phys. Rev. D 42 
(1990) 3321; M. F. Parry and A. T. Sornborger, Phys. 
Rev. D 60 (1999) 103504. 

[30] S. Kasuya and M. Kawasaki, Phys. Rev. D 61 (2000) 
041301; T. Hiramatsu, M. Kawasaki and F. Takahashi, 
JCAP 1006 (2010) 008; M. I. Tsumagari, Phys. Rev. D 
80 (2009) 085010; R. Allahverdi, A. Mazumdar and A. 
Ozpineci, Phys. Rev. D 65 (2002) 125003. 

[31] J. S. Liu, Monte Carlo Strategies in Scientific Comput- 
ing, Springer Publication (2001); C. P. Robert and G. 
Casella, Monte Carlo Statistical Methods, Springer Pub- 
lication (2010); M.Yu.Khlopov, Cosmoparticle physics, 
World Scientific(1999). 

[32] E. E. Flanagan, S. H. H. Tye and I. Wasserman, Phys. 
Rev. D 62 (2000) 024011. 



